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Abstract. A multipole spin-down equation based on a monopolar term is derived from the general expression 
i> = —f{i',t) and used to study pulsar evolution. We show that the time-independent version of such equation 
cannot reproduce the observed properties of pulsars and conclude that there is no equation of the form v = —fiy) 
consistent with the P — P diagram and braking index measurements. We explore the time-dependent model under 
the hypothesis of decaying magnetic fields, showing that an inverse linear decay gives reasonable evolutionary 
trajectories. This model distinguishes the evolution of Vela from that of the other three young pulsars considered. 
We discuss the origin of the monopolar term, which cannot be attributed to radiative processes, pointing to the 
importance of particle acceleration and/or mass loss processes in the dynamical evolution of pulsars. 
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1. Introduction 

Pulsars are accepted to be rapidly rotating and highly 
magnetized neutron stars powered by the loss of rota- 
tional kinetic energy. Yet, the mechanisms through which 
pulsars lose their rotational energy are not very well un- 
derstood. The discovery of pulsars was preluded by the 
rotating magnetic dipole model of Pacini I)19ti7|l . which 
predicted a spin-down relation v (y: , where v is the rota- 
tional frequency. Pulsar spin-down has then been studied 
through the power-law differential equation. 



(1) 



where fc is a constant and n is called the "braking in- 
dex" , determined by the physical mechanisms spinning- 
down the star. For a pure magnetic dipole n is equal to 
3 (Gold 119681 Pacini il968p . while a pure gravitational 
or electromagnetic quadrupole one has n = 5 (Ostriker 
& GunnHnSl Ferrari & Ruffini T^S^. The expansion 
of the Larmor equation into electromagnetic moments 
higher than dipolar gives n > 3. Departures from these 
basic models have been studied. For example, the de- 
formation of magnetic field lines by corotation with the 
magnetosphere might produce 1 < rt < 3 (Manchester & 
Tavlor [T977|l . pulsar winds n < 3 (Blandford & Romani 
119881 Manchester et al. I1985|l . while magnetic field de- 
cay (Chanmugam & Sang I1989|l and the alignment be- 
tween rotation and magnetic axis tend to give n > 3 
(Goldreich 1970). Physical processes such as phase transi- 



tions in the interior of the fast rotating neutron stars can 
produce changes in the moment of inertia / and led also 
to deviations from n = 3 even if the pulsar is spinning- 
down by pure dipolar radiation (Chubarian et al. I2U0UI 
Glendenning et al. I1997|l . However, these phase transi- 
tions are significant only for very fast pulsars, with pe- 
riods P < 1.5 ms. For pulsars with periods P 3 ms, 
these changes in the moment of inertia can be neglected. 
But while physical mechanisms to produce n < 3 have 
been proposed, the reconciliation of pulsar evolution with 
n < Z remains to be studied. This question, addressed 
to the case of young classical pulsars, is one of the main 
objectives of this work. 

The braking index n is frequently defined in terms of 
observational quantities, with a definition generalized to 
a second braking index, to: 
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Accurate measurement of both braking indices requires 
long-term timing measurements, a difficult task due to the 
existence of glitches and timing noise in young pulsars, af- 
fecting the measurement of i) and ly. Out of more than 
1300 detected pulsars, only a handfuU have reliable brak- 
ing index measurements (Table Still, measurements 
continue to be made and objects like PSR J1119-6127 
and PSR J1846-0258 might probably be recent additions 
to the list (Camilo et al. 2000. Mereghetti et al. 1^0(12)1 . The 
second braking index m is even more difficult to measure 
and has only been estimated for PSR B1509-58 (Kaspi et 
al. lBMI and the Crab pulsar (Lyne et al. .l988)l (Tabled). 
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Table 1. Timing parameters and stationary multipole 
model fitting parameters for the four young pulsars se- 
lected. 



Pulsar 


Crab'i) 


1509 - 58'2) 


0540 - 69'^' 


Vela''') 




OO.O 


1 KA n 


oU.o 


on o 






lOOO 


1 fifiA 


1 1 Qnn 
1 ioUU 


n 


2.509 


2.837 


2.01 


1.4 




±0.005 


±0.001 


±0.02 


±0.2 


m 


10.23 


14.5 








±0.03 


±3.6 






5(Hz-3) 


3.5 ■ 10'^"^ 


6.2 ■ 10'^^ 


6.6 ■ 10"!® 


3.8 • 10-1® 


r(Hz-i) 


1.0 • 10^1* 


1.6 ■ 10"" 


7.1 ■ 10"!^ 


1.3 • 10-1^ 


s(Hz) 


3.4 • 10-12 


2.0 ■ 10-12 


5.7- 10-12 


1.2 • 10-12 


Pbirth 


9.9 


14.2" 


17.6 


50 


{msr 










a p 

^birth 


defined as the period one 


dynamical time ago. 



^ Integration is stopped at tdyn = 1 year. (1) Lyne et al. 
rMEl : (2) Kaspi et al. (1994); (3) Manchester & Peterson 
JHHHJ; (4) Lyne et al. (il996l) . 



It is well-known that measured values of n and m are be- 
low, but close, to those of a magnetic dipole in vacuum 
(n = 3 and m = 15). 

Pulsar evolution is mostly studied through the P — 
P diagram, which shows a sparse distribution of pulsars 
clustered roughly at logP —0.3 and logP ~ —15. In 
principle it should be possible to link the bulk of the pul- 
sar population in the P — P diagram with the younger pul- 
sars through evolutionary tracks consistent with a general 
spin-down equation. However, neither the magnetic dipole 
model or a constant braking index equation can link these 
groups and it has been argued that pulsars might be born 
with small period derivatives, in a region of the P — P dia- 
gram where no pulsar detection has been reported (Camilo 
[l£96). Alternatively, magnetic dipole moments might de- 
cay during the pulsar life f Camilo 119961 Tauris & Konar 
IM?T)l . On these lines, Colpi et al. H2000|l studied the distri- 
bution and evolution in the P — P diagram of anomalous 
X-ray pulsars (AXPs), considering dipolar spin-down of 
magnetars. They conclude that magnetic field decay is re- 
quired to account for the observed distribution and X-ray 
luminosities. 

We propose here a spin-down model derived from the 
very general spin-down law v — —f{v,t). Although the 
derivation is empirical, the model accomodates the usual 
spin-down mechanisms, namely magnetic dipole radiation 
and gravitational radiation. In addition, a monopolar term 
is introduced, aiming to reconcile braking indices with 
trajectories in the P — P diagram. We calculate evolu- 
tionary trajectories for PSR B1509-58, PSR B0540-69, 
the Vela pulsar (PSR B0833-45) and the Crab pulsar 
(PSR B0531+21). The model is presented in Section 2, 
while in Section 3 we apply the model to the four pulsars 



just mentioned, assuming time-independent spin-down 
mechanisms. In Section 4 we study the model consider- 
ing time-dependent spin-down mechanisms. In Section 5, 
a discussion of possible physical mechanisms related to 
the spin-down terms is given, to lead to the concluding 
Section 6. 

2. The Model 

The basic assumption of this work is that the frequency 
derivative v can be described as a function of frequency v 
and time t only: 

y = -f{y,t). (3) 
We then assume the following properties: 

1. / > for all {u.t), as / describes energy and angular 
momentum losses. 

2. / is a continuous function of time. This requires that 
the cumulative effect of glitches on evolutionary paths 
in the P — P diagram can be neglected. 

3. / is an antisymmetric function, f{—v,t) = —f{i>,t). 
This is equivalent to ascribe a sign to the frequency v 
and its derivative i^: v and u have the same signs when 
the pulsar is spinning-up and opposite signs when the 
pulsar spins-down. This conditions will be further dis- 
cussed below. 

The model consists then of a simple Taylor expansion of 
/, restricted to the three lowest order terms: 

V = -~s{t)v - r{t)v^ - g{t)v^ . (4) 

Higher order terms can be neglected for the frequen- 
cies measured in known isolated pulsars i. A variable mo- 
ment of inertia, 1^0, would introduce an extra pos- 
itive monopole term I/2I that can be incorporated in 
s(t). This would affect only very fast pulsars and we will 
consider / constant in here. Equation |0J represents a 
very general spin-down model for isolated pulsars. Beyond 
its simple empirical derivation, it can represent the loss 
of rotational energy through physical processes described 
by standard electromagnetic and gravitational radiation 
multipolar terms, together with the ad- hoc "monopolar" 
term. We will argue later that this term has to be associ- 
ated with other type of processes, like particle acceleration 
or pulsar winds. We underline that if this model cannot 
be used to interpret the P — P diagram, then there is 
no analytical equation of the form given by equation 0, 
consistent with the required conditions, able to reconcile 
braking indices with pulsar evolution. 

3. The stationary multipole spin-down model 

Leaving momentarily aside the discussion of the origin of 
the different terms, we consider now the model given by 
equation when the coefficients g, r and s are constant, 

1 Binary and millisecond pulsars are not considered here. 
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following Alvarez & Carramifiana H1998fl . These coeffi- 
cients can be calculated if the frequency and both braking 
indices are known. Taking the first and second derivative 
of equation Q and using equations l|2Jl, one arrives to a 
matrix equation which can be inverted to give: 

/.9^VM)\ . / 3 -3 1\ / 1 \ 

rv^l{-y) = - -10 10 -2 n . (5) 

\svl{~v)j ^\ 15-7 \)\m-n^j 

The factors gv^/{—v), rv^/{—v) and sv/l—v) repre- 
sent the fractions of the total energy loss contained in 
the quadrupole, dipole and monopole terms respectively. 
Constraining each term to be between and 1 leads to 
the bounds on their values as function of n represented in 
Figure ni The value of m is restricted to the range 

Max(rt^ + 3n — 3; tl^ + 7n - 15) < m < + 5n - 5 , (6) 

while 1 < n < 5. The minimum and maximum allowable 
values for m are 1 and 45, which happen for n — 1 and n=5 
respectively. Accordingly, pulsars must be in a restricted 
region of the {n,m) plane defined by equation JHIl, and 
shown in Figure. |21 Spin-down evolution is initially domi- 
nated by the quadrupole term (n — > 5), afterwards by the 
dipole term (n 3) and finally by the monopole term 
(n — *■ 1). Under this model the braking index itself would 
be a rough qualitative age indicator. 

From the known timing parameters, v, i>, n and m, 
one determines {g, r, s} uniquely, from the matrix rela- 
tion ISJ, with no degree of freedom. This could only be 
done for PSR B1509-58, while for the Crab, PSR B0540- 
69 and Vela pulsars we assumed m within the constraints 
of the model. Evolutionary tracks were calculated for these 
pulsars, integrating equation Q backwards one dynami- 
cal time and forward to a point where the monopole term 
dominates the spin-down evolution (Fig. The results 
for each of the four pulsars can be summarized as follows: 

— the second braking index of the Crab pulsar, 
m = 10.23, is below the minimum value allowed by the 
model (to > 10.82) and in order to compute evolution- 
ary trajectories we used m — 11. This value of to is con- 
sistent with the initial period calculated in other mod- 
els for this pulsar (Mereghetti et al. .2002,. Glendenning 
I1996|) . The Crab is the only pulsar studied here with 
the precise age known, allowing some specific compar- 
isons. Integrating back to its actual birth (949 years 
ago) one finds an initial period ^18 ms (Fig0J. If we 
integrate back one dynamical time (1258 years) from 
now, we get an initial period of 10 ms, lower but within 
a factor of two of the period at birth. 

— for PSR B1509— 58 the measured values of to and 
n are consistent with the model. However, the inte- 
gration of equation Q could not be carried out one 
entire dynamical time backwards, as J dt converges 
at about 1350 years, i.e., before the present dynami- 
cal time (1553 years). Integrating back to a dynamical 
age of 1 year, we obtain a putative "initial period" of 
about 14 ms and a P two orders of magnitude higher 




1 2 3 4 5 



Braking index 

Fig. 1. L/Erot ratios for the quadrupole (upper panel), 
dipole [middle panel) and monopole [lower panel) terms. 
Only values inside the shaded areas are allowed by the 
stationary version of the model. When n —^ 1 and n — > 
5 the rotational losses are due only to the monopole 
or quadrupole term respectively. The observed braking 
indices give upper and lower limits for the different 
L/Erot ratios, as indicated by the dotted lines for Vela, 
PSR B0540-69, the Crab and PSR B1509-58 (from left to 
right in the horizontal axis). 

than for any of the other three pulsars. This pulsar 
seems to have been formed in a particular region of 
the P — P diagram. Its quadrupolar parameter g is 
at least two orders of magnitude higher that for the 
other pulsars. An initial phase of strong gravitational 
spin-down would naturally lead to the high P and P 
observed in this pulsar. 

— PSR B0540— 69: its second braking index is still un- 
known so we took TO — 8.0, in the middle of the range 
of allowed values. As the known timing parameters of 
PSR B0540-69 are similar to those of the Crab, albeit 
a somewhat longer period, their evolutionary paths are 
similar, with an initial P—P position of PSR B0540-69 
relatively close to that of the Crab. 

— Vela pulsar: the frequent glitches of Vela prevent any 
estimate of the second braking index, and render dif- 
ficult measuring its first braking index. We assumed 
TO — 3.5, consistent with the low first braking index, 
n = 1.4 ± 0.2 (Lyne et al. HM^ . As n < 2, Vela has 
already passed its minimum P value and its predicted 
evolution is close to a n = 1 line. We note that the 
track computed for the past evolution of Vela is close 
of that of PSR B0540-69. However, its motion in the 
P — P diagram is very slow and when integrating back 
one dynamical time Vela only reaches P « 50 ms with 
a braking index almost equal to 2. 
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Fig. 2. Allowed values of the braking indices n and to. 
According to the stationary version of the model, pulsars 
evolve from the upper-right to the lower- left, increasing 
their braking indices, but keeping them inside the triangle- 
like region. Both indices have been measured only for 
PSR B1509-58, which lies inside the allowed region, and 
the Crab pulsar, located slightly outside the allowed re- 
gion. For Vela and PSR B0540-69 the value of to has not 
been measured and we have assumed to to be in the mid- 
dle of the allowed range (open circles). Vertical dot lines 
indicate the value of n of these pulsars. Errors of n for the 
Crab and PSR B1509-58 are smaller than the dots. 




-3-2 10 1 

log,,, [Feriod/sJ 



Fig. 3. Evolutionary trajectories computed for the Crab, 
PSR B1509-58, PSR B0540-69 and Vela. The small dots 
on the tracks mark the present positions of these four 
pulsars, while the big dots mark the positions integrated 
back one dynamical time. The open circle marks the posi- 
tion of the Crab at the time of the SN1054 explosion. For 
PSR B1509-58 the integration converges before one dy- 
namical time, the large filled circle marking the P ~ P po- 
sition with tdyn = 1 year. The stationary model takes 
these pulsars into a region of the P — P diagram where 
pulsars have not been detected. 



Although exact tracks cannot be computed for older 
pulsars, their known timing parameters not constraining 
{g, r, s} sufficiently, we evolved backwards 17 pulsars with 
tdyn < 10^ years, assuming n — 1.5, in order to locate the 
region of their likely initial conditions in the P—P diagram 
(Fig. 13). We stopped the computation when tdyn = 1 year. 
We note that most of these pulsars arrive to higher P 
than all of the youngest pulsars, with the exception of 
PSR B1509-58. In order to get to P — P positions more 
consistent with that of the Crab, these pulsars would need 
to have a braking index closer to 1. A prediction of the 
stationary multipole model is that n <; 1 for the majority 
of the pulsar population. 

A more detailed discussion of the stationary model can 
be found in Alvarez (1998,) • Here we want to point out sev- 
eral inconsistencies of the model, which lead us to reject 
its non-time dependent version. As it can be seen in Fig.|31 
the evolutionary tracks of the youngest pulsars take them 
to a region of the P — P diagram where no pulsars have 
been found. Although a large P might difficult the de- 
tection of pulsars located there, the region corresponds 
to pulsars slower but more energetic than some known 
radio-pulsars, specifically those with same P but lower P. 
Pulsars in that region should be detectable. We believe 



these tracks are indicative of the inconsistency between 
the model and the observations. A second discrepancy is 
the eventual dominance of the monopole term, which im- 
plies vjv constant (FigE)). As n — > 1 the dynamical 
time tends to a constant, t^yn — > l/2s, which values be- 
tween 2.8 X 10^ years for PSR B0540-69 and 13 x 10^ years 
for Vela. The dynamical ages of the rest of the pulsars in 
the P — P diagram are much larger and cannot be repro- 
duced by the model from the data of these young pulsars, a 
clear inconsistency between the predicted evolution of the 
young pulsars and the timing parameters of older pulsars. 

Through the study of the stationary multipole model 
proposed here we reach the general conclusion that there 
is no time-independent evolutionary equation v — —f{y) 
able to reconcile braking indices with the P — P diagram. 
Note that including & v on term'^ cannot prevent the 
eventual dominance of the monopolar term, needed to 
explain the braking index of Vela (and probably also 
of PSR B0540~69). The impossibihty to have a time- 
independent expression means that the physical proper- 
ties involved in the dynamics of pulsars must vary on 
timescales comparable to their life as active radio-pulsars. 
The most obvious form of time dependence is through the 

^ inconsistent with the antisymmetric condition. 
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Fig. 4. Pulse period as function of time predicted for past 
and near-future evolution of the Crab pulsar. The evo- 
lutionary trajectories are given by the all allowed values 
of m. The dashed, full and dot-dashed lines correspond 
to m =13.64, m =11.0 and m =10.82 respectively. The 
vertical dotted lines mark the present time {t = 0), the 
moment of the SN1054 explosion {t = —948 yrs) and one 
dynamical age backwards. 

magnetic moment of the star, namely magnetic field decay 
or alignment between the magnetic and rotational axes. 
We consider in the next section the time-dependent mul- 
tipole model using physically expected time-dependences. 

4. The time-dependent multipole model 

Modelling the distribution observed in the P ~ P di- 
agram requires considering plausible time dependences 
in physical parameters involved like the magnetic dipole 
moment of the neutron star. Observational evidence in 
favor of magnetic field decay is discussed by Tauris & 
Konar (^2001,) , who analyzed the relation between char- 
acteristic time against rotational energy of the pulsars 
in the Princeton Catalog (Taylor et al. I1998|l . Additional 
evidence comes from the measurements of their proper 
motions and from the study of ultramagnetized neutron 
stars or magnetars. Sang & Chanmugam fTUUU) argued 
that dynamical and kinetic ages of pulsars can only be 
reconciliated with decaying magnetic fields of the form 
B{t) = B{0)/{l+t/tc)- Further evidence of magnetic field 
decay is provided by anomalous X-ray pulsars (AXPs) in- 
terpreted as magnetars. Colpi et al. H2UUU|I required the 
magnetic field in these neutron star related objects to de- 
cay as a power law on timescales of the order of 10"* years. 

We introduce time dependence in our model through 
the coefficients {r{t),s{t)} oc B{Q)'^4'{t/'tc), where tc is 
the characteristic time-scale for field decay and ip a di- 



Fig. 5. Evolutionary trajectories for 17 middle aged pul- 
sars evolved backwards to tdyn = 1 year. The four 
youngest pulsars (Crab, PSR B1509-58, PSR B0540-69 
and Vela) are denoted by the large dots. 

mensionless function satisfying "0(0) = 1. Frequently used 
functional forms of ip are exponential or inverse linear. For 
simplicity, we will only consider future evolution for the 
youngest pulsars, so that the quadrupolar coefhcient g(t) 
can be neglected hereafter^. We can then study evolution- 
ary tracks on the P — P diagram writing equation in 
terms of the period and its derivative as 

P=(^+,soP) ^{t/Q. (7) 

This equation describes the period evolution of pulsars 
for a decaying magnetic field. It can be integrated analyt- 
ically, leading to: 

P^(t) = {ro/so + Pq) exp [2sotc "^{t/t,)] - ro/so , (8) 

where Pq = P(0) and ^-(2;) = J^^ ip{u)du. If -0 -> for 
t ^ tc,'^ will generally converge to a value '^00 and the pe- 
riod will tend to a constant, P^o- This occurs because the 
spin-down mechanisms disappear as "0 approaches zero. If 
one knows the present values of ip and 0, the initial pa- 
rameters ro and sq can be derived from equation Q and 
its time derivative: 

where n is the present braking index and a decay-law 
means tp < 0. The decay timescale tc is taken here as a free 

^ Neglecting g also means that we do not require knowing 
m. 
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5 10 15 

t (1000 years) 

Fig. 6. Dynamical versus chronological time for the Crab 
pulsar. Note that tdyn 4.6 x 10'^ yrs for ages greater 
than ^ 20 X 10'^ years. As the model forces O (x ly, as 
n — > 1, dynamical ages tend to constant values which, for 
the four pulsars studied, are below the values observed for 
other pulsars. The dotted lines mark the date of birth of 
the pulsar and today. 

parameter to fit the data, to be determined by selecting 
physically meaningfuU trajectories in the P — P diagram 
and To > and Sq > 0. Note that rg and sq are both 
positive if 



l>i^>n-3, 
P/P 



(10) 



which relates tc with the dynamical time and requires 
n > 1 always. On the other hand, old pulsars can now 
have arbitrarily large braking indices, as n can go to oo 
when tjj/tp ^ cte, P Poo and P ^ 0. The largest loss of 
rotational energy occurs when ip —i- 0, ^ constant i.e. 
when n is close to one. 

In the next sections we present evolutionary tracks for 
the cases of exponential and linear magnetic field decay. 

4.1. Exponential decay 

Based on the observed deficit of pulsars with periods larger 
than 1 second, Ostriker & Gunn H1969|l proposed a rapid 
decrease of pulsar radio luminosity as result of exponen- 
tially decaying magnetic fields. 



B{t) ^ Bo exp i-t/tc) 



(11) 



a natural result of ohmic dissipation. In this case the 
model becomes 



Table 2. Time-dependent model parameters and expo- 
nential magnetic field decay timescales consistent with 
tq > and sq > 0. 





Pulsar 




ro 




So 






ic(yrs) 












) 


(Hz) 




min 


best 


max 




Crab 




1.0- 10" 


14 


4.0 - 10" 


V2 


3300 


4.5 ■ 10" 


8.0 ■ 10" 




1509 - 


58 


2.1 ■ 10" 


13 


1.8 - 10" 


12 


3390 


3.7 ■ 10" 


1.5 ■ 10^ 




0540 - 


69 


1.0 - 10" 


14 


6.8-10" 


12 


6600 


2.5 • 10" 


5.0 • 10" 




Vela 




2.3-10" 


17 


1.7-10" 


12 


1.1 • lO'^ 




2.0 • 10^ 



PSR B0a33-'I5 



PSR B1509-58 




P : 



-2t/ta 



(ro/P + soP) 



(12) 



-1.5 -1 -0.5 0.5 1-1.5 -1 -0.5 0.5 1 
I'Og jj h'orlod] Log jjHcrlod] 

Fig. 7. Evolutionary trajectories using exponentially de- 
caying magnetic fields for the Crab, PSR B1509-58, 
PSR B0540-69 and Vela. The three curves correspond to 
the lowest, best and largest tc consistent with observa- 
tions. For Vela the best trajectory is that of minimum 
tc, as shown in Table 13 Trajectories were followed for 
lOtc, etc and 5tc for the PSR B053H-21, 5.5*^ and Ate 
for PSR B0833-45, lOOtc, 5.5tc and Ate for PSR B1509-58 
and lOOic, 6*. and A.8tc for PSR B0540-69. 



and we have 'foe — 1/2. For the youngest pulsars we ob- 
tain the trajectories shown in Figured Satisfactory tra- 
jectories are obtained for tc ~ 10'* yrs (Table EJ, which 
correspond to extremely rapid evolution, as the magnetic 
field becomes negligible in just 10 tc and pulsars cross 
the death-line of the P — P diagram in lO** to 10^ years. 
This timescales are too short to be consistent with esti- 
mated birth rates of one pulsar every 90 years (Brazier 
& Johnston I1999|l and supernova rates of one every 10- 
30 years (Van der Bergh & Tammann I1991|l . We con- 
clude that exponentially decaying magnetic fields cannot 
be used to properly fit the P — P diagram without con- 
tradicting present known pulsar and supernova rates. 
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Table 3. Time-dependent model parameters and inverse 
linear magnetic field decay timescales consistent with vq > 
and So > 0. 



Pulsar 


ro So 
(Hz-i) (Hz) 


ic(yrs) 
min best max 


Crab 
1509 - 58 
0540 - 69 
Vela 


9.4 • 10"'=' 7.5 • 10"'^ 

1.2 • 10~" 4.8 • 10~^^ 

7.3 • 10"^=^ 1.0 • 10"" 
2.1 • 10"^* 1.7 • 10"^^ 


2090 1.0 ■ 10* 3.0 ■ 10" 
1830 1.0 ■ 10* 4.0 ■ 10* 
4980 1.0 ■ 10* 2.0 ■ 10* 
1.0-10^ 



4.2. Inverse linear decay 

mm 



Sang & Chanmugam 1)1990(1 proposed an inverse linear 
decay law for the magnetic field of pulsars, 



B{t) = 



BjO) 

1 + t/tc 



(13) 



to reconcile dynamical and kinetic ages. Using this 
decay law in the model, relatively short timescales, 
tc ~ 10"* years, are needed to fit the evolutionary trajec- 
tories of the four pulsars considered with the P — P distri- 
bution, as shown in Figure|H|for the Crab, PSR B1509~58, 
PSR B0540-69 and Vela. Except for the Vela pulsar, to 
be discussed below, trajectories can be fitted through the 
bulk of the pulsar population. Although the timescales 
needed to fit the data are similar to those for the case of 
exponential decay, linear decay is far more convenient as 
pulsars remain active for 10^ — 10^ tc, crossing the death 
line of radio emission in some 10^ years, in reasonable 
agreement with pulsar and supernova rates. 

For the Vela pulsar, the minimum decay timescales 
consistent with the model, i.e. with positive ro and sq, 
do not manage to turn its evolutionary tracks towards 
the bulk of pulsars. However, the uncertainty in the mea- 
surement of n allows us to consider a range of values. 
Evolutionary trajectories with n < 1.4 take Vela to a re- 
gion of large periods and high period derivatives where 
no radio pulsars have been found. For n = 1.2, Vela takes 
about 2.0 X 10^ yrs to slow down to P ~ 10 s, and 3.6x10^ 
years for n = 1.3. Taking these numbers and considering 
one pulsar like Vela born every 11000 years, we should 
expect to see between 18 and 27 pulsars in the same 
region of the P — P diagram occupied by Vela. Indices 
1.5 < n < 1.6 can produce reasonable tracks (Fig.|5Jl, lead- 
ing to longer periods than those of the other three young 
pulsars, but consistent with long-period pulsars present in 
the P — P diagram. It does not seem feasible to reconcile 
Vela trajectories with those of the Crab, PSR B1509-58 
and PSR B0540-69, which pass though the bulk of the pul- 
sar population without impossing a braking index n > 1.6, 
inconsistent (at least at the la level) with observations. 
So, even though the model can fit the data, the evidence 
points to the Vela pulsar as a particular type of pulsar. 
Pulsar evolution might distinguish between Crab-like pul- 
sars, with magnetic-fields decaying at periods Poo ^ 1 s. 
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Fig. 8. Evolutionary trajectories for the Crab, 
PSR B1509-58, PSR B0540-69 and Vela, using an 
inverse linear decay-law for the magnetic field. The 
decay timescale used is tc — 10'' except for Vela where 
tc = 1.04 X 10^ years. The trajectories were run for 10'^ tc, 
except that of Vela which was run for 10<c- 



and -less common- Vela-like pulsars, where the magnetic- 
field decay is slower and becomes spun-down to asymp- 
totic periods between 1 and 10 seconds. 
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Fig. 9. Evolutionary trajectories for the Vela pulsar, con- 
sidering different possible values oin. Trajectories using n 
larger than ~ 1.5 are consistent with known long-period 
pulsars, while n < 1.4 is inconsistent with observations. 
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5. Physical Spin-down Mechanisms 

The main feature of the model presented here is the in- 
clusion of the monopolar term, required to derive reason- 
able evolutionary tracks considering the measured brak- 
ing indices. The other two terms of the model have been 
widely used in the past and are usually interpreted in 
terms of electromagnetic or gravitational radiation losses. 
An inclined rotating magnetic dipole in vacuum produces 
a dipolar term given by 



E, 



dip 



3c3 



(14) 



where VL ~ 2t:i' is the angular frequency and ii± the 
component of magnetic dipole moment perpendicular to 
r2. Pacini (1967) was the first to equate Edip = 
introducing the basic equation of pulsar dynamics. A 
quadrupolar term can be associated either to the next 
electromagnetic moments or to the lowest gravitational 
radiation multipole moment: 



E, 



quad 



32G 



(15) 



where le the product of the moment of inertia and equa- 
torial ellipticity of the star. The fact that observed indices 
are close to 3 indicates that gravitational radiation is not 
the dominant form of energy loss in pulsars at present, 
although it might be important in the very early stages of 
pulsar evolution. 

Directly equating the Larmor equation to rotational 
losses cannot lead to the monopolar term proposed in 
equation Q . The Larmor equation arises from considering 
the energy flux of an electromagnetic wave, S oc \Erad\'^k, 
with the radiation field Erad proportional to the accelera- 
tion of the charged particles, scaling therefore as fi" with 
n even and > 4. Even in the framework of the models men- 
tioned in Section 2, it is unlikely that radiative emission 
processes can give a different dependence on 17. 

On the other hand, there is clear evidence that par- 
ticle acceleration and massive winds are important en- 
ergy loss processes in pulsars. This evidence includes the 
short lifetimes of relativistic electrons in the Crab Nebula 
and the need of a source to account for its luminosity 
(as pointed out first by Oort and Walraven ll956|l . obser- 
vational changes directly observed in the vicinity of the 
Crab pulsar (Hester et al. I1995|l and the X-ray images of 
the Crab, PSR B1509-58 and the Vela pulsar (Weisskopf 
et al. HUM Gaensler et al. "MI^, and Helfand et al. IMITI 
respectively). The monopolar term cannot be of radiative 
origin, so we propose here its relation with particle and/or 
mass loss processes. We note we cannot discard a. E oc il'^ 
term but this term would be neither necessary nor suffi- 
cient to explain the data; therefore we avoid introducing 
it and assumed the antisymmetric condition for simplic- 
ity. Neglecting the E oc fl'^ term (equivalent to i> cx j/^) 
we put forward the issue of whether particle acceleration 
and/or mass-loss follow an antisymmetric spin-down equa- 
tion, like the Larmor equation does. 



A monopolar spin-down term arises from any energy 
loss mechanism with a. E oc fl^ dependence, which can be 
expressed as 



En 



(16) 



where /3 is a dimensionless factor that can depend on time 
but not explicitly on f2. It is an open question whether 
energy-loss processes like particle acceleration and pulsar 
winds can provide this sort of expression and we can only 
especulate about this point. Michel l,1969)i estimated the 
torque exerted on a neutron star by a pulsar wind as 



T = (7r$V4c) 



(17) 



where $ is the magnetic fiux carried away by the wind. 
Assuming $ can be constant, a pulsar wind would lead 
to a monopolar term in the energy loss equation. A de- 
tailed analysis of pulsar wind energetics has been given 
by Harding et al. |1999), who assumed an electric current 
flowing in the magnetosphere causing a magnetic torque 
on the star. The corresponding loss of rotational energy 
due to this torque is given by: 



TP — 1 n2 7-)6 



1 



/?2 /?2 



(18) 



where Rlc — c/U,. Ropen is the distance from the center 
of the star to the magnetosphere region at which the mag- 
netic fleld lines become open due to the pulsar wind. For 
Ropen independent of fl, this implies E^ind cx 17^, i.e., a 
monopolar-likc term. 

6. Conclusions 

A multipole equation for pulsar evolution can be derived 
naturally as the Taylor expansion of the general equa- 
tion z> = —f{i',t). While it might give insight into short 
term evolution of pulsars, the time independent version of 
the multipole equation is inconsistent with observations, 
indicating that no equation of the form i> = —f{v) is 
able to fully model pulsar spin-down evolution. A time 
dependent multipole equation considering exponentially 
decaying magnetic-flelds cannot reproduce the data with- 
out shortening the active life of pulsars below 10^ years. 
Introducing an inverse linear magnetic field decay with 
timescales of the order of 10'* years gives good agree- 
ment with observations and allows pulsars to remain ac- 
tive about 10^ years, in agreement with pulsar and super- 
nova birth rates. The model is consistent with observations 
of the Vela pulsar only if its braking index is ^ 1.5. The 
different trajectories of this pulsar relative to those of the 
Crab, PSR B1509-58 and PSR B0540-69 suggest the exis- 
tence two types of pulsars: those with relatively short de- 
cay timescales which evolve to periods Poo ^ 1 s and Vela- 
like pulsars, with more persistent magnetic fields evolving 
to the larger period pulsars observed in the P~P diagram. 
We believe the need to introduce a monopolar term in the 
spin-down equation is an indirect observational evidence 
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that particle acceleration and/or mass- loss processes must 
follow EocQ'^. 
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